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Astronomy, 



The Casimir force and free energy at low temperatures have been the subject of focus for some 
time. We calculate the temperature correction to the Casimir-Lifshitz free energy between two 
parallel plates made of dielectric material possessing a constant conductivity at low temperatures, 
described through a Drude-type dielectric function. For the transverse magnetic (TM) mode such a 
calculation is new. A further calculation for the case of the TE mode is thereafter presented which 
extends and generalizes previous work for metals. A numerical study is undertaken to verify the 
correctness of the analytic results. 



PACS numbers: 72.20.-i,11.10.Wx,42.50.Lc,78.20.Ci 

There has been an explosion of interest in the Casimir 
effect [H, generalized to dielectrics by Lifshitz [2], since 
the modern experiments began with Lamoreaux in 1997 
Q. The zero-temperature Casimir-Lifshitz theory seems 
to have been confirmed to 1% accuracy over a range from 
100 nm to a micrometer [1, i, i, 0, S S [M [H [ll • 

However, there has been a continuing controversy over 
the temperature dependence of this effect. The pre- 
scription given in Ref. was seriously questioned by 
Bostrbm and Sernelius [l4| who pointed out that neces- 
sarily the transverse electric reflection coefficient at zero 
frequency must vanish for metals. This discontinuity pre- 
dicted a new linear temperature term at low tempera- 
tures, resulting in about a 15% correction to the result 
found by Lamoreaux. Lamoreaux believes that his ex- 
periment could not be in error to this extent p^ . More 
heatedly, Mostepanenko and collaborators have insisted 
that this behavior is inconsistent with thermodynamics 
(the Nernst heat theorem), because it would predict, for 
an ideal metal, that the free energy has a linear tempera- 
ture term at low temperature, and hence that the entropy 
would not vanish at zero temperature |l6j |. Moreover, 
they assert that the precision Purdue experiments rule 
out the linear temperature term in the low temperature 
expansion 12J. 

The issue is as y et unresolved, and is summarized in 
recent reviews jlTl . ITsj . We will not add further to the 
discussion of this controversy here. Rather, the purpose 
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of this paper is to examine another purported tempera- 
ture anomaly. In several recent papers [l^ [U, HI] 
Geyer, Klimchitskaya, and Mostepanenko have claimed 
that in real dielectrics, which possess a very small, but 
nonzero conductivity which vanishes at T = 0, a similar 
discontinuity in the transverse magnetic reflection coef- 
ficient occurs, which would lead to a similar violation of 
the Nernst theorem. The same applies to semiconduc- 
tors whose conductivity vanishes as temperature drops 
to zero. The solution according to these authors, as in 
the TE case for good conductors, is to prescribe the ef- 
fect away. We argue, however, that such a solution is 
physically unsatisfactory. 

In Sec. [H we will review and clarify their argument for 
a standard Drude-type permittivity model for a weakly 
conducting material. We thereafter work out the leading- 
order temperature corrections to the free energy in the 
cases where the media are assumed to have a finite but 
small residual conductivity at T = 0, as is implied when 
a Drude model is employed for taking the conductivity 
into account. This is a new result to our knowledge 
(a similar calculation for materials with zero conduc- 
tivity was undertaken in Ref. [l^). While this calcu- 
lation does not solve the thermal anomaly brought forth 
in Refs. [H, [13, HH and reviewed in section |T1 it serves 
to further illuminate the mathematical behaviour of the 
free energy of poor conductors at very low temperatures 
when different models for the dielectric response of the 
materials are employed. A similar calculation is subse- 
quently performed for the TE mode, which extends that 
of Ref. 2^ in several ways: We allow for the conduc- 
tivity to be small; we work out one further order of the 
temperature correction to the free energy; and we allow, 
for generality, the permittivity to have a finite dielectric 
constant term in addition to the Drude-type dielectric 
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response due to free charges. 

A word about units. For our theoretical calculations, it 
is most convenient to use Gaussian electromagnetic units, 
as well as natural space-time units: ?i = c = fce = 1- 
However, for final results, which could be experimentally 
observed, we use SI units. The mapping between units 
is very simply carried out by dimensional considerations, 
using the unit conversion factor he — 1.97 x 10~^ eV cm. 
The conductivity transformation between Gaussian and 
SI units is the simple replacement 47r<T = cr^^/eo, where 
Eq = 8.85 • 10~^^ F/m is the absolute permittivity of the 
vacuum and the notation used to explicate that SI 

units are used. 



I. TEMPERATURE ANOMALY FOR 
SEMICONDUCTORS 

Here is a simple way to understand the argument of 
Ref. [l^ . Suppose we model a dielectric with some small 
conductivity by the permittivity function 
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The essential point is that asC— >0,£^eif(T = 0, 
otherwise e oo. The Casimir (Lifshitz) free energy be- 
tween two halfspaces, separated by a distance a, assumed 
to be of the same material for simplicity, is given by 
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F = -y 

m=0 



ln(l - ' 



(1-2) 



Here K,rTE, and ttm are functions of the discrete Mat- 
subara frequencies Cm — 2TrmT; = fc^-|-C^ with kj_ the 
transverse wave vector, directed parallel to the surfaces. 
As is conventional, the prime on the summation mark 
implies the m = term be taken with half weight. We 
need to examine the behavior of the reflection coefficients 
in the small C limit. These are 



r-TE 



rTM 



K+ C^(£ - 1)' 



(1.3) 
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where e = e(iC)- For the case of an ideal metal, it was 
which was discontinuous: 



rTE(C-0) = 0, 



lim rxE = -1, 



(1.5) 



so this gave a linear temperature term when the sum 
over Matsubara frequencies is converted to an integral 
accor ding to the Euler-Maclaurin formula, for example 
(Ref. '19| uses the Abel-Plana formula, but that is equiv- 
alent) . 



For a dielectric the TE reflection coefficient is con- 
tinuous and vanishes as C ~* Oj but if there is a small 
(but not zero) conductivity which vanishes with T lin- 
early or faster, the TM coefficient exhibits a discontinu- 
ity at C = as we now explain. When the conductivity is 
small we can assume there exists a temperature so that 
the m = 1 Matsubara frequency, Ci — ^ttT, satisfies the 
inequality 



< Ana < Ci < Wo 



in which case 

rTuiiC = 0) 



1, rTM«l) 



e+ 1 



(1.6) 



(1.7) 



Typical values of luq are in the optical or near IR fre- 
quency regions, so Eq. (jl.6p will hold at room tempera- 
ture for many semiconductors. If now a goes to zero as 
T ~> linearly or faster, Eq. (|1.6p continues to hold true 
all the way to zero temperature where it becomes a true 
discontinuity: 



rTuiiC = 0) = 1, 



lim rTM«) 
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(1.8) 



Clearly if a reaches some residual value > 0, Eq. (|1.6p will 
not hold near zero temperature. Likewise the discontinu- 
ity disappears should a be exactly zero in a temperature 
region of finite width including T — 0. 

As in the metal case, Eq. (|1.8p gives rise to a linear 
temperature term in the pressure and the free energy 
(see e.g. Ref. [iTj and references therein for details). Let 
fm be the summand of Eq. (|1.2p or a similar expression 
for the Casimir pressure. Since fm is discontinuous at 
TO = 0, we must replace it by a continuous function. 
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where f^n is continuous 

fo- 



fni = IhUm-o/m; TO = 
fm — fm 1 ^ 7^ 



(1.10) 



so that the Euler-Maclaurin summation formula can be 
applied to the sum over /„. Then the first term in the 
third form in Eq. Ijl.9p gives rise to a free-energy con- 
tribution which is a linear function of T . Defining the 
shorthand notation 
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that hnear term is 



F™ = — / dK K[ln(l - Aoe-^-^") - ln(l - e'^-^")] 
47r Jo 
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[Li3(Ao)-C(3)] 



(1.12) 
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where the polylogarithmic function is 



fc=i 



(1.13) 



Note that the hnear term vanishes for e ^ cxd as is clear 
from noting the relation to the Riemann zeta function: 

Li„(l) = C(")- (1.14) 
Thus at zero temperature, the entropy is nonzero, 

1 



dT 167ra2 



[Li3(Ao)-C(3)], (1.15) 



which, if physical, is a violation of the Nernst heat the- 
orem, or the third law of thermodynamics, which states 
that the entropy of a system must vanish at zero temper- 
ature. 



II. GENERAL FORMALISM 

For reference throughout the next sections we will go 
though the formalism of determining the leading temper- 



ature corrections to the Casimir (Lifshitz) free energy by 
use of the Euler-Maclaurin formula, a procedure often 
employed previously. 

Considering one polarization mode at a time, the free 
energy for the q mode {q =TM,TE) is written in the form 



= /(a,T)E g{m) 



(2.1) 



m=0 



where we have pulled out a convenient prefactor. 

When T — + the Matsubara sum becomes an integral, 
so the temperature correction to the free energy, given 

by 



AF, = /(a,r) 



m=0 



dm 



g{m), 



(2.2) 



can be determined by use of the Euler-Maclaurin formula. 
For the summands of the Lifshitz formula, the higher 
derivatives of g{m) are singular near m — 0. When this 
is the case the Euler-Maclaurin formula can be applied 
to the sum starting ai m — 1 (or a higher value of m) 
instead, whereby 



r = 



E 



dm 



g{m) 



i.9(0) - g{m)dm + ig(l) - £ (Syi^^^'^-^' (1 



^5(0) - 9im)dm + ^g{l) ~ ^g'{l) + - 



(2.3) 



where Bn are the Bernoulli numbers. 



are now 



B2 = -; ^4 = — B(i = —-... (2.4) 
6' 30' ^42' ^ ' 



using the convention of [2J| §23.2. [Two remarks are 
called for here: We have assumed that g and all its deriva- 
tives vanish at infinity, and we have converted this for- 
mula into one which is commonly asymptotic because we 
have omitted the remainder term which is present when 
only a finite number of derivatives terms are retained. 
Thus we are considering only the leading terms in an 
asymptotic expansion for small T.] 

As mentioned above, g{m) is not analytic at to = 0. It 
can be written in the asymptotic form for small to 

5 (to) -Co + cim + cg/aTO^^^ 

-|- C2;TO^ InTO, + C2TO^ + . . . , m ^ 0. (2.5) 



The terms needed for the right hand side of Eq. 



5(0) - Co, (2.6a) 

5(1) = Co + Ci+C2 + Cy2 + (2.6b) 

g'il) ^ ci + C2i + 2c2 + -C3/2 + . . . , (2.6c) 



g"'{l) = 2c2, - -C3/2 + . . . , 



dmg{m) 



1 1 1 

Co + - gC2/ -I- -C2 



rC3/2 



(2.6d) 



(2.6e) 



When inserted into Eq. (|2.3p the terms involving cq and 
C2 cancel and one is left with 



~ ci , lies; 49 

1 W Cq/o 

12 360 1920 ' 



(2.7) 



Here the term due to Ci is exact, whereas the terms refer- 
ring to logarithms and half-integer powers of to, C21 and 
C3/21 receive contributions from all higher derivatives in 
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the Euler-Maclauriii formula, and to obtain exact expres- 
sions for the coefficients, ah such terms must be kept, as 
we now show. 

Retaining the higher derivative terms in the Euler- 
Maclaurin formula one finds by using 



^3/2 



m — l 



3(4n - 7)! 
"24n-5(2n-4)!' 



and 

i>2n 



J2n-1 



-TO In TO 



2(2n-4)!, 



m— 1 



dvn?"- ^ 

that the temperature correction to free energy is 

AF, =/(a,T)f, 
where with Eq. ((O)) 

*C2Z + $C3/2 + . . . , 



n > 2, 
(2.8a) 

n > 2, 
(2.8b) 
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with the coefficients 



B2 
2 

2 

5 ' 



\ ^ B2n1p2n 
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E 

n=2 



(2«)! 



(2.9) 
(2.10) 

(2.11a) 
(2.11b) 



These series are formally divergent as is typical for per- 
turbation series near singularities. Indeed, they arise 
from the asymptotic Euler-Maclaurin formula (|2.3p . For 
example, $ can be recognized as a special case of the ex- 
pansion of Riemann zeta function in terms of Bernoulli 
numbers, Eq. 23.2.3 of Ref. [2^ with an infinite number 
of terms retained in the sum, and the remainder omitted. 
A meaningful value can nonetheless be assigned to them 
through Borel summation as detailed in Appendix [X] 
Numerically, the mathematical software Maple computes 
the numerical values by means of a Levin u-transform to 



* = 0.03044845705840. 
$ = -0.0254852018898. 



(2.12a) 
(2.12b) 



By either numerical or analytical correspondence we thus 
recognize that 



$ = C(-|). 

C(3) 
47r2 ' 



(2.13a) 
(2.13b) 



where C is the Riemann zeta function. 

When Eq. (|l.ip is used in the Lifshitz formalism with 
constant and finite a and £ > 1 in Sees. IIIII and IIVI 
we will find that the terms of F stemming from ci, C3/2 
and C21 are proportional to , T^/^ and respectively. 
Higher-order terms of 5 (to) will likewise give higher-order 
temperature corrections. 



III. TM MODE, RESIDUAL CONDUCTIVITY 

In the following sections we will work out the low tem- 
perature behaviour of corrections to the free energy un- 
der the assumption that a Drude-type dielectric func- 
tion (jl.ip may be used, and that a is finite and constant 
with respect to ( and T for small T and As argued 
in Ref. [25|, when a is finite close to zero temperature 
Nernst's theorem will be satisfied. Here we will calcu- 
late explicitly the low temperature behaviour of the free 
energy for the TM mode. 

Conventionally, semiconductors are found within 
the broad interval of conductivity a in SI units 
10"5(r2m)-i < crSi ^ io5(f]ni)-\ that is 



10^s~^ < CT^Veo < 10 



16„ 



(3.1) 



For numerical purposes we will use the intermediate 
value fj/eo — lO^^s^^, which is large enough not to 
hamper numerical verification unnecessarily, but small 
enough to distinguish the material in question from a 
good metal. The frequency corresponding to a/eo for a 
metal is ujp/v, where ujp is the plasma frequency and ly 
the relaxation frequency. For gold at room temperature 
i^p/iy has the approximate value 3.5 • lO^^s^^. 

Returning to Gaussian units, we consider the TM mode 
and introduce the shorthand notation 



t = 



and the symbol 



2ttT 



Ana Ana h{a^^/eQ) 



fx — mt. 

If Ana = 10^2 s""'^ as assumed above, 
t « 0.83T 



(3.2) 



(3.3) 



(3.4) 



with T in Kelvin. 

The free energy is given by Eq. p.2p , for which we now 
consider only the TM term: 
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F™^^y / d«;«;ln(l-^e-2-), (3.5) 



m=0 

where the reflection coefficient squared is 



^ — ' TM — 



(3.6) 



Here and henceforth the index to on Matsubara frequen- 
cies Cm and quantities dependent on it will frequently be 
suppressed. 

The temperature corrections to the free energy at low 
temperatures are dominated by small frequencies, so we 
can assume as an approximation that the middle term of 
Eq. (jl.ip is simply equal to e — 1 and write 



e«) K£ + 



Ana 



1 



(3.7) 
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We define the dimensionless quantity 



2a 



a = 2a(47rcr) = — {a^^/eo), 



(3.8) 



where a is the distance between the semiconductor plates. 
For the value Ana ~ lO^^s^"'^ or smaller, a is a small quan- 
tity, which we use to define a criterion for the smallness 
of the conductivity in the remainder of this paper: 



a < 1. 



(3.9) 



For a — Ijim and a as above, as used for numerical pur- 
poses later, a has a value of about 6.7 • 10^^, so this 
criterion is well satisfied. 
By defining the variable x 



X = 2Ka = 



na nafj. 



Ana 



c 



(3.10) 



A can be written 



A 



1 + £^ - ^yr+[r+(£^^n3/i]a2^ 



(3.11) 

and the integral p.Sp with the use of Eq. (|3.2p and ( — 
2'KmT becomes 



F 



TM 



_ {AT:aft 



m=0 



where 



/>oo 

g{m) = / dxa;ln(l - Ae~^). 



(3.12) 



(3.13) 



We wish now to extract explicitly the temperature de- 
pendence of the integrals g{m) in Eq. (|3.13p . The proce- 
dure we choose is to expand Eq. p.l3p to leading order 
in the small parameter a, and then expand the resulting 
term in powers of m to get the form (|2.5p . 

The first term in the Taylor expansion of the logarithm 
in powers of a is 



ln(l - Ae-^) = -Lii(A^e-^) -f 0{a^) 



(3.14) 



where we use the polylogarithmic function defined in 
Eq. (|1.13p and define the quantity 



A,, 



l-f(£-l)/i 

l + (e + l)/i 



(3.15) 



For integral s < 1 the polylogarithm Lis(j/) can be ex- 
pressed by elementary functions, specifically 



Lii(y) = -ln(l - y); Lio(y) 



Li-i = 



y 



(3.16) 



The summand ^(m) thus has the form 



g{m) 



f 

J (XL 



dxx\Xx{A^e^'')^0{o^). (3.17) 



Now we will expand gim) in powers of m. It is easy to 
show from Eq. (|1.13p that 

dy\Jx.n{Ce-^y) = -iLi„+i(Ce-'^^); (3.18) 

from which by partial integration 

g{m) = -a//Li2(A^e-"'^) - 03(^^6""^) + O^a"). 

(3.19) 

We now use the property 

Li„(Ce-y) = Sr^Li„_KC) (3.20) 



for \C\ < 1 to expand the polylogarithms in powers of 
afx. The terms containing Li2 then cancel and we are 
left with 

gim) = -Li3(A^) + ^a^ti^Lh{A^) + ©(a^) (3.21) 



with Afj_ given by Eq. (|3.15p . Henceforth we shall denote 
the first two terms of the expansion (j3.2ip giijn) and 
gii[m). The remaining 0{a^) term comes from the error 
in Eq. (|3.14p . As before we are going to truncate the 
expansion in a at leading order, but will evaluate the 
explicit correction ~ to (|3.2ip later as a measure of 
the error. We thus have the simple expression 



gi{m) = -Li3(A^). 



(3.22) 



We will next expand Eq. (|3.22p in Li3(A^) does not 
have a Taylor expansion near m = (where Aq — 1) 
because its second derivative is singular here. Using 



d 1 
— Li„(?/) -Li„_i(y), 
dy y 



we differentiate Eq. (|3.22[) to find 
4tLi2(A^) 



(3.23) 



(3.24) 



[l + {e + l)^i][l + {e~l)^i\ 
We can use the identity |i26|] 



Li2(z) + Li2(l-z) = — -ln(z)ln(l-z), (3.25) 
6 



which is easily verified by differentiation, use of Eq. p.l3p 
and Li2(l) = 7r^/6. Furthermore Li2(l — A^) has a simple 
Taylor expansion around = 1, 



and 



Li2(l-Ap) =.4/i-4(e + l)^2^ 



4 - 8e/i 



[H-(£--f l)/i][l + (e-l)M] 



(3.26) 
(3.27) 
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whereby we find 



a'lim) = 



^ 4mtM ^+4) +16m^2ln4/i- 



(3.28) 

where the next term of the series is of order t^. 

Comparing with Eq. (|2.5p we recognize the coefficients 



2nH 



Cl 



C21 = 8t' 



which we insert into Eq. (|2.10p to find 



E 

.m=0 



dm 
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(3.29) 



Sm^. (3.30) 



We thus obtain the approximate correction to the free 
energy for small t: 



^ A 2 2 



E 



dm 



gi{m) 



(47rcr)3 



f [tt^ - 144*t] 



(3.31) 



in terms of our reduced units t and a. In SI units insert- 
ing (I2l3bl) : 



TT^ikjiT)^ ^ C(3)(fcBT)3 



72?i(aSVeo)a2 7r[;i(aSi/eo)a]2 

nHkBT)' ( _ 72C(3)fcBr ' 
'72?i(aSVeo)a2 1^ ^^ha^^/eo 



(3.32) 



A. Correction due to subleading terms of 
Eq. ([33T|) 

Twice in the above we truncated the expressions at 
leading order in the parameter a, in Eq. p.l4p and 
Eq. (|3.2ip . As an indication of the magnitude of the 
error we will calculate the next order in a of Eq. p.2ip 
while a similar calculation for Eq. p.l4p is more trou- 
blesome due to singularities and beyond the scope of the 
present effort. The correction cx of Eq. (|3.2ip was 

Ag(m) = ^a^m^t^Ui{A^) + 0{a^t^). (3.33) 

We will only consider the first term, since the next terms 
give temperature corrections cx and higher. We Tay- 
lor expand as before in powers of /i 

Lii(A^) = -ln(4t) -Inm-I- (e + 2)/i + e'(Ai2), (3.34) 

wherewith the leading correction from Agi{m) is found 
from Eq. (j2.10p to order to which only the term cx 
him contributes: 



(3.35) 



Being a independent, the correction (|3.35p is much 
smaller than the leading term (|3.3ip for small a. In SI 
units: 



SF 



TM 



C(3)(fcBT)^ 



(3.36) 



The relative magnitude of this term compared to the 
term of Eq. p.32p is with our numerical data 



2„2 



(^"Veo) ^a 
4c2 



2.8-10" 



(3.37) 



The correction from the truncation of Eq. ()3.14p is likely 
to be of similar size and therefore much smaller than the 
accuracy of the numerical investigation. 



B. Numerical investigation of TM mode result 

The numerical investigation in Fig. [1^ employs 
Eq. dni) with cr^Veo = lO^^g-i^ ^ ^ ^1.67 and = 
8 • lO^^s"^ as appropriate for Si [27]. While the analyt- 
ical expression fits well for T < O.IK, corrections ~ 
become important beyond this point. The two leading 
orders in temperature corrections were shown to be inde- 
pendent of e to leading order in a, but the correction 
(not calculated analytically herein) depends heavily on 
this value. A qualitative measure of this effect is given 
in Fig. [TJ) where we have used £ = 1, cetera paribus. In 
all plots the theoretical curve is that given in Eq. (|3.3ip . 
ignoring corrections. 

It is noteworthy that, as seen from Fig. [U while the 
term of Eq. (|l.ip gives the dominant temperature 
correction for small T, nearly all (99.7% with our data) 
of the free energy at T = is due to the e term. 

While the fit pictured in Fig. [1] is indicative, a much 
more sensitive confirmation of the accuracy of the theo- 
retical results is provided by considering the quantity 



R = 



AFth - AF„ 



AFt 



th 



(3.38) 



where AFnum is the direct numerical calculation and 
AFth is the theoretical result to next-to-leading order, 
in the form p. 311) . An analysis exactly like this was per- 
formed in Ref. |23>] ; the reader may refer to that paper 
for further details. 

We have found that AFth is of the form 



AFth = -CT^l - CiT) 
and assume AFn„ni to be of the form 



AF„„„, - -DT\l - DiT + D2T' 



from which one finds 
C -D D 



(3.39) 



(3.40) 



R 



C 



(Cl - D^)T 



C 

-^[D2 + Ci{Ci-Di)]T^ 



(3.41) 
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FIG. 1: F"^^ and its approximation with (a) e = 11.67 and 
(b) e = 1. Correction curves, calculated Irom Eq. (|3.3Hl . are 
shifted to match the numerical calculations at T = in each 
graph. 
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FIG. 2: The quantity R defined in Eq. ((338|l plotted for the 
TM result (|3.3H) and numerical calculations. 



Fig. [5] to Eq. (I3.42P it is clear that D2 < which implies 
that the coefficient of the term of the free energy be 
positive, which conforms with the corrections in Fig. [T] 
not accounted for to order T^. 



IV. TE MODE, RESIDUAL CONDUCTIVITY 

For the TE mode the dominant temperature correction 
to the free energy comes from the last term of Eq. (jl.ip . 
The permittivity (jl.ip . which can be approximated as 
Eq. (|3.7p . is similar, but not identical, to that for a Drude 
metal, considered in Refs. [H, [l^. There, instead of 
Eq. (|3.7p the permittivity was assumed to be 



UJZ 



In the special case where C 
becomes 



D and Ci — Di, this 



R = -D2T^ + 0{T^) 



(3.42) 



which is zero in the limit T — Q and has zero slope in 
this limit. The zero temperature limit of R and its slope 
thus provide measures of the accuracy of the theoretical 
results: if the coefficient is correct, R should approach 
zero as r ^ 0, and if the coefficient is correct, the 
slope of R{T) should vanish in this limit as well. We 
have not taken the corrections l|3.36p into account in the 
plotting of Fig. [51 

We have undertaken a numerical study of the be- 
haviour close to zero temperature, resulting in the graph 
of R shown in Fig. O Due to the vanishing denomi- 
nator of Eq. (|3.38p . the analysis is extremely sensitive 
to numerical errors as the zero temperature limit is ap- 
proached. From the figure it seems clear that the errors 
in the two coefficients are small enough to confirm the 
correctness of Eq. (|3.3ip , although some caution must be 
exerted due to the numerical volatility of R. Comparing 



^-metal 



C(C + ^) 



1 + ^ 



E. 



(4.1) 



The principal difference is that the constant term e is 
assumed to be significant here and kept general. Since 
for small C the term ^ dominates the constant term, 
an approximation to the low-temperature behaviour of 
the dielectric would be expected to be found by the same 
analysis as that of Refs. [2J,[2^ but with the substitution 



— Ana. 

V 



(4.2) 



For typical semiconductors, 47r(T is smaller than io^jv for 
a good metal by many orders of magnitude. For this rea- 
son, since the free energy at zero temperature is of the 
same order of magnitude for the metals and semiconduc- 
tors for the same separation, the relative temperature 
corrections for the TE mode are expected to be much 
smaller than for a metal. Thus, there is reason to inves- 
tigate whether the effects of e > 1, while negligible for a 
metal, could be important for small tr. In some dielec- 
tric materials, as is well known, e can exceed unity by 
as much as two orders of magnitude, and a more careful 



8 



analysis is therefore justified. The procedure is the same 
as above, and an extension of that found in Ref. [23|], to 
which the reader may turn for further detail. 

It was found in Refs. (H,!!! that for T ^ 0, and e = 1, 



TE 



where 



C,= (i^(21n2-1), 



C5 ^ 



48 
2^ 



6 



C(-3/2)(47ra)3/2a. 



(4.3) 

(4.4a) 
(4.4b) 



Here (^{y) is the Riemann zeta function (for this closed 
form of Cs , see Appendix A of Ref. [23]). 

For the numerical values indicated this gives the SI 
values 



C2 = 1.6185719 • 10 



Cs = 2.5844373 • 10 



-19 



Km^ V 1012s 



(4.5a) 



-22 



f a \ f cr^Veo 



K3m2 VlQi^s-i 



(4.5b) 



Thus the TE temperature correction is expected to be 
positive and in the order of magnitude of lO^^^J/m^ at 
T = IK. 

The numerical calculations shown in Fig. [3] were com- 
plicated by the fact that the thermal corrections are 
many orders of magnitude smaller than the free energy at 
zero temperature, making a graph of the quantity R sim- 
ilar to Fig. [2] unfeasible within the assumption of a ^ 1. 
We show here, however, that assuming £ > 1 does not 
change the theoretically predicted thermal correction to 
the free energy to order T^, and therefore merely refer to 
Ref. 23] for further numerical support of the theoretical 
result. 



A. More general treatment of the TE mode 

Let us treat the TE mode temperature correction to 
the free energy more carefully. Starting with the expres- 
sion (|1.2p we perform the substitution 



K KjJL 

" Cv/^(iCF^ ~ xC 

where we define the recurring quantity 

Then the free energy may be written 
TE _ (47ra)3t . " ' 



F 



47r2 



■JI 5(w) 



(4.6) 



(4.7) 



(4.8a) 
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. 5 

Temperature (K) 



FIG. 3: Temperature dependence of the free energy for TE 
at 1 micron separation. The solid line is an exact numerical 
calculation including all terms of Eq. with e = 11.66 and 
u>o = 8.0T0^^s~^, the dashed line is the parabolic temperature 
correction (|4.3p . The term oc is too small to be visible 
in the graph. 



with 



g{m)=x I dxxln{l- Be-"^''). 



(4.8b) 



The squared reflection coefficient given by (II. 3p now de- 
pends only on x: 



B = {x 



\fxp- 



We expand the integrand of g{m): 



ln(l - 56""^"=) = ln(l - B 



axxB 
l-B 



(4.9) 



(4.10) 



Note that this is as far as we can expand this way, 
since the next term of the a expansion gives a divergent 
contribution (an alternative method which avoids some 
divergences but is somewhat more cumbersome is the 
method employed in Appendix A of Ref. [23j where the 
corrections are calculated without the use of the Euler- 
Maclaurin formula.). 

Consider the first terms of the expansion ()4.10p 
(we dub the terms of g{m) from the expansion 
9i{m),giiim), ...): 



gj{m)=X^ dxx\n[l- {x- \hF+Ty]. (4.11) 



This integral can be evaluated explicitly (a similar in- 
tegral was evaluated in Ref. [231] where the lower limit 
was approximated as zero). Perform the substitution 
X — sinhw. Then we may write 



gi{m) 



XI 

4 



du(e""-e-^")ln(l-e-*") (4.12a) 



9 



with 



X 10 



uq — arsmn— = — in — , — 

2 \^Fim-^ 



X 

With the substitution y 

.2 rvo 



(4.12b) 



-2 m 



gi{m) 



X 



dy{y-'-l)\n{l-y'), (4.13a) 







where 



2/0 



^JepL + 1 + ^ 
= l-2V7^+2^+(e-2);u3/2 



(4.13b) 



The integral is straightforward to evaluate and the result 
is 

,2 



91 (m) 



X_ 

8 



1 

2/0 



+ yo ln(l - yl) 



-2yo + 21n 



1 



2/0 



1 - 2/0 



(4.14) 



We expand this in powers of /i and find that the terms 
oc /^'^/^ cancel, consistent with the small-x dependence of 
the integrand of gj{m). We are left with 

2 

gj{m) ^(21n2- 1) - ^[ln4^ + e(21n2 - 1)] 



Comparing with Eq. (j2.5p we see 
ci = —-(2 In 2 — 1) and C21 



(4.15) 



(4.16) 



while the dependence on e only enters in the C2 term 
~ m? which does not contribute to the Euler-Maclaurin 
formula. The temperature correction to first order in a 
is thus 



47r2 



2 In 2 - 1 
48 



m 
T 



(4.17) 



We see that the leading term conforms with Eq. (|4.3[) 
when Eq. (|4.4ap is inserted. The first term beyond those 
calculated is proportional to T^/^ according to (|4.15p . In 
SI units with (|2.13bp : 

AFTB^^!!(W!(,i^2-l)-M%^. (4.18) 



48eo?ic2 



The term of (14.18^ has the same form as that found 
for ideal metals in the limit aT <C 1 [29l [30| . A similar 
term is present in Eq. (|3.36p . [Note that the correc- 
tion for the TM mode is not fully accounted for therein] . 
A numerical comparison of this term with the difference 
between the graphs in figure [3] is shown in figure 2) It 
shows that the coefficient in (|4.18p is the right order of 
magnitude, but the numerical precision is not sufficient 
to draw definite conclusions about its accuracy at this 
time. 



15 - 
10 - 



F - F 

quad num 

?(3) (k T) V (Snh^c^) 



0.1 0.2 0.3 0.4 0.5 

T (K) 



FIG. 4: The difference between the numerically calculated 
free energy and the quadratic T term of Eq. (|4.18p (equal to 
the difference between the graphs in figure [3)l plotted against 
the absolute value of the term of (I4.18|l . 



B. First order correction to expansion (|4.10p 



The first order correction term in Eq. (|4.10p is easily 
calculated with a similar scheme. We have 



1~B 



gii{m) =ax^ / dx 



ax^ ^^ x{x-Vx^f 

4 J^/^ V^^^TT 



(4.19) 



The procedure for solving this integral is as before. Sub- 
stitute X = sinh u to get with a little shuffling 



(4.20) 



With the substitution z — e ^ this becomes very simple: 



with 



(4.21) 



zo = e 



/ V£At + 1 - Vm 

V y/sn + 1 + ^ 



1/2 



Thus we find the /i expansion of gii[m): 
5//(™) = f 0/i^/^-(2-£V/^ + . 
Hence, with Eqs. ((2T0l) and ((48aP . 



TE _ (47rcr)^a 5/2 



II 



487r2 



$^5/2 _^ 0(^7/2)^ 



(4.22) 



(4.23) 



(4.24a) 
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or in SI units with $ = C(— 1): 



-|)a(aSVeo)3/2 



6h 



3/2 



5/2 



(4.24b) 

Comparison with Eq. (j4.4bp shows fuU agreement with 
the result for metals (e — 1). 

It is worth noting that while the next-to-leading tem- 
perature correction is of order T^^'^, the term cx T"^ in 
Eq. (|4.18[) dominates it with respect to a. Thus in the 
small a limit the T^/^ dependency becomes all but im- 
perceptible. 



V. CONCLUSIONS 

We have worked out the two leading terms of the tem- 
perature correction to the Casimir-Lifshitz free energy 
at low temperatures between poor conductors obeying a 
Drude-type dispersion relation. We have assumed that 
the materials have a small residual conductivity (com- 
pared to the reciprocal of the interplate separation) which 
is finite and constant with respect to temperature and 
frequency near T = 0. 

The calculation for the TM mode complements that 
of Ref. [l^ where the temperature correction for free en- 
ergy between two dielectrics of zero conductivity was cal- 
culated. Both the TE and TM free energy temperature 
corrections are quadratic to leading order. To the extent 
of our computations, the TM mode has integer powers of 
T beyond the leading whereas the TE mode has both in- 
teger and half-integer powers. The temperature anomaly 
reviewed in Sec. [J occurs when the transition from finite 
to zero conductivity happens at exactly T = 0, and while 
the analysis presented here does not resolve the anomaly, 
it is of interest to know the low temperature behaviour 
of the free energy in each of the two cases. 

Note furthermore that the effects of the static dielec- 
tric permittivity e only enters to order for the TM 
mode and order T''^^ for the TE mode. The fact that 
the coefficient of the term T^/^ appears to depend on e is 
noteworthy since only integer powers of T were reported 
in Ref. [10|, although seeing as we have not calculated 
the coefficient here it is possible that cancellations occur. 

Our calculations are delicate since they rely on the 
relative smallness of different parameters simultaneously. 
We have assumed the parameter t (essentially temper- 
ature T divided by conductivity a) small while at the 
same time letting a be small compared to the inverse of 
the separation a. This is the reason why the leading or- 
der temperature corrections in (|3.3ip appear to diverge as 
a vanishes. On a deeper level these subtleties stem from 
non-commuting limits in the Lifshitz formalism which are 
the cause of anomalies such as that reviewed in section H] 
Another curious property both of the present calculations 
and those of Geyer, Klimchitskaya, and Mostepanenko 
figj is that the free energy corrections of order in T just 
beyond what we have considered here appear to diverge 



as e ^ 00, as indicated for example by Eqs. (|3.34p and 
(|4.23p . This limit would a priori be expected to yield the 
ideal metal limit. Such phenomena should be addressed 
in future studies in the effort to achieve full understand- 
ing of the low temperature behaviour of the Casmir force 
and free energy. 

The asymptotics of the Lifshitz formula as frequency 
and temperature approach zero are fraught with inher- 
ent subtleties both mathematical and physical. While the 
method employed herein is highly useful for its simplic- 
ity and transparency, it has limitations because the func- 
tions involved are not analytic in the limits considered 
and non-integer powers and logarithms enter. Physically 
we have assumed herein a model which may represent 
certain physical systems, but avoids the temperature be- 
haviour which leads to the anomaly reviewed in Sec. ID It 
also neglects effects which may be of importance, such as 
spatial dispersion, a subject which has been extensively 
investigated over the years [Sl, Ji3, 34]. A theoretical 
effort to attempt to describe the screening effects and di- 
electric response of the vanishing density of free charges 
in insulators near zero temperature involving all impor- 
tant physical effects will likely be required in the future 
and will hopefully provide the resolution of the anomaly 
for dielectrics. 
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APPENDIX A: BOREL SUMMATION 

The Borel sum of the divergent series X^J^o exists 
(Ref. H, Sec. 8.2) if the function 



E 



(Al) 



is convergent for sufficiently small x and the integral 

B{x) = / dte-*<j){xt) (A2) 
Jq 

exists. Then the Borel sum is J2'^=o ~ ^i^)- Consider 
the quantity ^' defined in Eq. (|2.11ap and consider the 
term 



B2„(2n-4)! ^Bn{n-4) 



(2n)\ 



(A3) 
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(the latter equality follows from i?2n+i = 
1,2,...). Letting a„_4 ~ Bn{n — 4)!/n! we get 



0, n 



n=0 ' 



1 



where the identity J2^=o BnX^/nl = xj {e^ 
The Borel sum 



X X 

2 ~ 12 

(A4) 
1) was used. 



^ =B{1) = 



dt( 



t4 



t 



t__t^ 

T " ^ 2 ^ 12 



(A5) 



is now possible to evaluate analytically (the divergence 
in the lower limit is illusory because the expression in 
brackets is of order t'^) to find the desired value as ^ = 
1/36 — 2^P — C(3)/(47r^). A similar numerical procedure 
using (|Aip and (|A2|1 will give the value of <i>, which as 
noted in the text has a well-known asymptotic expression. 
An alternative and equivalent approach which is most 



often simpler is to sum each term of the expansion of 
g(m) in (|2.5p directly (disregarding the zero temperature 
contribution) and obtain finite values of the terms of the 
temperature expansion by means of zeta regularisation 
(36| . Such a procedure immediately yields $ — C(~i) 
by the definition of the zeta function as the analytic con- 
tinuation of C(s) = X]m=i Likewise the value of ^ 
can easily be found by comparison with the asymptotic 
series expansion of the derivative of C,{s): 



In 



n. 



(A6) 



n=l 



whereby ^ = —(,'{— 2) — C(3)/47r^. The same reasoning 
also yields the coefficient of the ci term of (|2.10p directly 
as C(~l) = —1/12. We thank Emilio Elizalde for alerting 
us to this point. 
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